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If} . Abstract 
CN 

We study the finite-size scaling of heavy-light mesons in the static limit. We compute two- 
point functions of chiral current densities as well as pseudoscalar densities in the e-regime of 
heavy meson Chiral Perturbation Theory (HMChPT). As expected, finite volume dependence 
turns out to be significant in this regime and can be predicted in the effective theory in terms 
of the infinite-volume low-energy couplings. These results might be relevant for extraction of 



heavy- meson properties from lattice simulations. 



1. Introduction 



The simulations of heavy-light mesons made out of a heavy quark (charm or bottom) and a 
light one (up, down or strange) on the lattice are challenging because they require very large 
volumes in order to keep systematic errors under control. The reason is that the dynamics 
of these systems involve very distinct energy scales: the heavy-light (hi) meson mass, Mhi, 
the light pion masses M\\ and Aqqd, that should all be kept sufficiently below the UV cutoff 
(i.e. the inverse lattice spacing), and sufficiently above the infrared one (i.e. the lattice box 
size). Both requirements can only be met in very large lattices. 

If the heavy quark mass is sufficiently large a good effective description is provided by 
heavy quark effective theory (HQET) [H [21 [3] , which is obtained in the limit of infinite heavy 
quark mass, or static limit. In this limit, the scale Mhi disappears from the problem and the 
UV cutoff can in principle be as low as the cutoff used to describe light meson dynamics. 
Indeed this approximation has been extensively used to simulate heavy-light mesons in lattice 
QCD (for a recent review on heavy flavour phenomenology from lattice QCD see [4]). 

Whether the heavy quark is treated in the static limit or not, an obvious question is if we 
can do better concerning the constraint on the box-size. After all, the finite-size scaling of 
heavy-light systems should be dominated by light pions physics, since these are the lightest 
modes in QCD. To the extent that pion physics can be described by chiral perturbation theory 
(ChPT) , it is conceivable that finite-size scaling of heavy-light systems can be accurately 
predicted using ChPT, as the finite-size scaling of light mesons is O [7]. 

In this paper, we investigate the possibility to predict the finite-size scaling of heavy-light 
systems, when the lightest pions are light compared to the inverse box size, from chiral 
perturbation theory. We will consider this problem in two limiting situations depending on 
the mass of the heavy quark: 

• The heavy quark is significantly above the light one, but still treatable in ChPT: this 
would correspond to considering hi mesons in the the mixed-regime introduced in [8]. 

• The heavy quark is static and therefore chiral dynamics can be treated in Heavy Meson 
Chiral Perturbation Theory (HMChPT): this would correspond to considering hi static 
mesons in the e-regime. 

Even though these two situations are physically very different, the pion dynamics responsible 
for the finite-size scaling properties should be pretty much the same. It is therefore interesting 
to see explicitly how a quantitative matching of the finite-size effects takes place, by comparing 
the finite volume dependence of correlation functions in ChPT and HMChPT. 

We consider the two-point function of left-handed current densities that will be computed 
to next-to-leading order in the e-regime in both effective theories. We will also consider the 
two-point correlator of pseudoscalar densities to the leading order, since finite-size effects are 
important already at this order. Anticipating the possible use of these results in simulations 
we also present the results in the partially-quenched (PQ) case. 
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The paper is organised as follows. In section [2] we present the results for the two-point 
functions in the mixed-regime of ChPT, when the heavy quark is treated in the p-regime and 
the light ones in the e-regime, that is in the so called mixed-regime. In section [3] we discuss 
the formulation of HMChPT in the e-regime and present the results for the same correlators. 
In section [7] we compare both results and discuss the implications. In section [8] we briefly 
comment on the applications to lattice QCD and conclude in section EJ 



2. Heavy-light mesons in the mixed-regime of ChPT 

The goal of this section is to study the finite-size scaling of heavy- light mesons in ChPT, 
when the light quarks are in the e-regime. We assume that the meson is composed of a heavy 
quark of mass rrih and a light quark of mass mi and that both masses are very different 
mi <C rrih, but both can still be treated in the context of ChPT, that is 

Ml y = {mx + ™^ «(4nFf, x,y = h,l. (2.1) 

Under this hypothesis the finite-size effects at NLO are predictable by using the common 
ChPT Lagrangian, that is: 



F 2 

-ChPT — 



dJJdprf - ?Tr M ] U + rfM 
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(2.2) 



plus the counterterms one has to consider at one loop that were found by Gasser and 
Leutwyler and are proportional to the Low Energy Couplings (LECs) Li.[9\. The pseudo 
Nambu-Goldstone bosons are parametrised by U G SU(N), with N = Ni + Nh, being Ni 
(Nh) the number of light (heavy) quarks. We have absorbed the vacuum angle 9 in the light 
quark masses. That is, the mass matrix M is: 



ie 



M = {m x e N i , . . . ,m Nl e N i ,m Nl+ i, . . .,m N } . (2.3) 

The mesons are placed in a box of volume V = L 3 T, which is sufficiently large to contain the 
typical QCD scale, and the heavy meson mass scale, but small compared to the lightest pion 
mass: 

M hi L » 1, M U L < 1. (2.4) 

In this situation it is expected that the finite volume effects associated to the scale Mhi are 
exponentially suppressed, while those associated to Mn are not. This regime of ChPT has 
been named mixed-regime in [HI [TO], since some of the quarks are in the p-regime and some 
in the e-regime. A convenient power-counting for the quark mass and momentum in this 
situation is 

m/~ e 4 , m h ~ e 2 , L' 1 ~ p ~ e , (2.5) 
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so that the LO mass of the heavy-light mesons will be: 



Ml = ^ ~ e 2 

h — p2 



(2.6) 



We refer to [H [10] for further details on the implementation of ChPT in the mixed-regime, 
both in the full and partially-quenched theories. We just remind here that in this regime it 
is convenient to parametrise the pion field U like: 



U 



U 
1 



'111 

e f , 



with the perturbative pion field £ satisfying the condition: 

J d 4 xTr [T a Z(x)] =0 if T a e SU(Ni), a = 1, ■ ■ ■ , N? - 1. 



(2.7) 



(2i 



In this way, the light zero modes to be treated non perturbatively are collected in Uq, and 
consequently they are dropped from £. In the references [SI [10] two different parametrisations 
where used. We have explicitly tested the results with both of them. 

In refs. [SJ [TU], the light-light meson correlators were computed. Here we extend this 
computation to the two-point correlation functions of heavy-light left-handed currents and 
pseudoscalar densities, to relative 0(e 2 ) order: 



Tr [T a T b ]Cj{t) 



d 3 x(jg(x)J b (0) 



(2.9) 
(2.10) 



Tr [T a T b ]C P (t) = Jd 3 x(p a {x)P b {0)y 

where in QCD the current and pseudoscalar densities can be formally defined as 

J^ = ^T a llM P^, P a = i^T a l5 i,, (2.11) 

t represents the Euclidean time, x = (x, t), and P_ = (1 — 7s)/2. In order to represent a 
heavy-light meson, T a is any traceless generator with one index in the light subsector and 
the other one in the heavy one, for example: 



(T a ) 



;(5ih8ji + SuSjh)- 



(2.12) 



As usual, in ChPT these operators can be represented to leading order in the momentum 
expansion by 



F 2 

J" = — Tr 
n 2 



T a Ud^ 



P a = i-Tr 
2 



T a (U-U^ 



(2.13) 



1 In order to simplify notation, we use throughout the same notation for the operators in QCD and their 
representation in ChPT. 
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These results are useful in their own right to describe for example kaon correlators in a finite 
volume, when the s quark is in the p-regime and the u and d are in the e-regime. We will also 
be interested in isolating the finite volume effects that survive in the static limit rrih — > oo, 
which should match those obtained in HMChPT. In order to recover the results for various 
full and partially-quenched situations of interest we consider the following computations: 



Degenerate heavy quarks: the Goldstone manifold is SU(Nf l + iVj), with JV/, quarks 
of mass rrih and JVj quarks with masses mi (I = 1, • • • ,N{), with the counting rules 
of Eq. (|2.5p . We can consider then the quenched limit of the heavy quarks by taking 
the replica limit Nh — * 0. These results should match in the rrih — * oo limit those of 
HMChPT, where the heavy quarks are treated as static sources and all the light quarks 
are in the e-regime. 



Non-degenerate heavy quarks: the Goldstone manifold is SU(Nh + N s + iVj), with 
Nh — ► quarks of mass rrih (i.e. the valence heavy quark), N s sea quarks of mass m s 
and N of masses mj 4 , where both rrih ~ m s ~ e 2 . This can be matched to HMChPT 
in the limit rrih oo. This situation corresponds to having sea quarks both in the e 
and in the p-regimes, for example if one considers B or D mesons in 2+1 dynamical 
simulations, where the s quark is in the p-regime and the u and d quarks are in the 
e-regime. We can also quench the light quarks N — > (quenching the heavy sea quarks 
N s — > is equivalent to Case A) , which would then correspond to the study of D and B 
mesons in a PQ mixed-action approach with sea quarks in the p-regime and the valence 
light quark in the e-regime. 

2.1. Conventions 

We describe in the following our conventions for the propagators that we use to write down 
the results in a compact form. 

The propagator for a pion with mass M is, in finite volume: 



Since some zero modes are factorised in the mixed regime, we also need to consider propaga- 
tors in which they have been subtracted: 



• Case A. 



• Case B. 




(2.14) 




(2.15) 
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The singlet part of the propagator gives rise to the following functions: 

(r> 2 ) 2 F(r> 



E(x,N s ,N h M) _ v g (p 2)2 F ( P)iVs)iVz)M ) NfVM 2 ' 



with □ = d Xfl d Xfl , and 



1 e *p^ 

UE(x,N s ,N h M) = > - ON -, (2.16) 

1 j F^(p 2 )F(p,7V s ,A^M) 1 ; 



N Ni 

F( P ,N s ,N l ,M) = ^L^ + ^. (2.17) 



Once we integrate over space, the correlators exhibit exponential decay at large distances. 
This is represented by the function: 

[ 1 , s COSh \ M (I" - 1*1)1 

F( '- M)S / dxG( *' M)= 2Msl\4] ' <2 ' 18) 
when the pion running in the line has a mass of order e 2 , or by: 

T fcl (l)=/*xB(,,0)-l[(M_>) , _^ 
if the mass is of order e 4 . 

When two mesons propagate we need to introduce the function: 



(2.19) 



1 f HP HP 

k 00 (M 1 ,M 2 ,t) = -J2<2—(t,M lp ) — (t,M 2p ) 

p I 

- (p(t,M lp )^-{t,M 2p ) + (Mi <- M 2 )j 1, 



where we have introduced the shorthand M ap = ^/ M 2 + p 2 . This expression is substituted 
by: 

fcooCMx,*) = A lim o ^fcoo(Mi, M 2 , t) + ) , (2-20) 
when M2 lays in the e-regime. 

2.2. Left-current correlator 



Case ^4 
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The result for the left-correlator at NLO using the above definitions is: 
F 2 

C { j A \t) = ^M? A) P(t,M {A) ) 



+ (Ni-^\k O0 (M h ,t)\ 1 (2.21) 



where we have defined 



Ff A) = F 2 - I [N h - ±j G (0, M hh ) - (0, 0) - ^±^G (0, M h ) 

' + ^A G ^ M vh) + lE(0,N h ,Ni,M hh ) + SM^{2L A N h + L B ), (2.22) 



2(N l +N h ) 2N h J v " hJ 2 



Mf A) = M 2 h 



1/9, s 2N h + 3N, , 

1 - ^ ( %M 2 h {2UN h + L 5 - 42V A 2* - 2L 8 ) - * ' 2 G(0, M,J 



! Dg(M^V^l4^ \ 2/1 iVj t A 

+ 6 J " Tk \m ~ T " 4 ((C/o + ^MJ 



(2.23) 



and 



Hi = rmZV, (2.24) 



M„ 2 , = —J—M* (2.25) 

while M^ is defined in Eq. f|2.6j) . 

A few observations are in order. The UV divergences in i* 1 ^ and in M 2 A ^ can be shown 
to cancel in the renormalisation of the NLO couplings of Gasser and Leutwyler, L\s. We 
have also checked that the result matches the result of [8] for non-degenerate quarks in the 
e-regime in the appropriate limit. 

This result represents the finite-size scaling of kaon-like states (m^ = m s and mi = m u = 
rrid) in the mixed-regime for various situations: 

• 2 + 1 dynamical simulations setting: = 1, Ni = 2, 

• PQ simulations where the h quarks are quenched and the I quarks are dynamical by 
taking the replica limit N h -> of Eq. (f2T2Tj) . 

• PQ simulations where the I quarks are all quenched or partially quenched, while the h 
quarks are dynamical. In this case, the appropriate value of Ni must be taken, but also 
the zero-mode integrals ((Uq + Uq)u) need to be properly defined d. 



2 Note that one cannot consider a fully quenched theory with Nh = Ni = on the basis of Eq. (I2.21[) . 
because the singlet has been integrated out [11] . 
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We discuss now the result of the zero-modes integrals ((Uq + Uq)u) (for further details see 
[12\ I10j). In order to treat the situation where some light quarks might be quenched, we 
distinguish within the light (e-regime) sector Ni sea quarks and N q quenched ones. When 
restricting to a topological sector v, the averages in all the cases described above can be 
obtained in a compact and general form, in terms of the partition functional |13^ 114]: 

crv n \\ - det [/ x j~ 1 *^+j-l(^)]ij=l,-2JV g +JV 1 ,„ 

*N q ,Ni+NMm) - N q , 2 2 2N q+Nl , 2 2 . ' (2.2b) 
Llj>i>l\f 1 j H'i) LLj>i>N q +l\t J 'j H'i) 



and its derivatives. Here ^7's are defined as Jv+j-i^i) = (—iy K v +j-.i(fii) for i = 1, • • • Nq 
and J u+ j-i([jLi) = I u+ j^\(fj,i) for i = N q + 1, • • • 2N q + N\, where I u and K v are the modified 
Bessel functions. For the observable of interest here, the result in the theory with all the 
light quarks dynamical is: 

{{Uo+ 2 U ^ = Aln^ +0 ( {W} ), (2.27) 
while in theories where the valence light quark is quenched is: 

^ + ? )w) " - lim ^yln^^XM). (2-28) 



Case B 

In this case we will denote respectively the squared mass of the heavy-light mesons and the 
decay constant at NLO by M? B -, and F(b) ■ The result in case B in the replica limit N^ — > 
is: 

C ( J B \t) = ^Mf B) P(t,M (B) ) 

T [ ( 1 \ — M 4 N N 

-2V\ [ Nl -Wj^^ + Nl{ Ml h N- 



MLN + M* h N - 2Ml h MlN 



k 00 (M h , M hh , t) + N s k 00 (M s ,M hs ,t) 



{Ml h N-MlNf 

_MUMl h -Ml) fj koo(MhiM2it) ) 1 (2 . 29) 



M 2 =M hh 

where we have defined the shorthand 

r2 _ N l 



M n = Ar ' (2.30) 
Vs N+N s ss v ; 
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and N 



N s + N h while 



(B) 



F 2 - ^G (0, M hs ) - ^G (0, M.) - §G(0, 0) - ^G(0, M h ) + ±^(0, N„ N h M ss ) 




Ml 



(B) 



We have performed several consistency checks of these results. For m s = mh Case A is 
recovered. UV divergences do cancel. We recall that G(0, 0) has no divergences in dimensional 
regularisation, and it is given by |15j : 



where (3\ is a so-called shape coefficient, which depends on T/L. 

In the replica limit iVj — > 0, this result represents the finite-size scaling of kaon-like correlators 
in PQ simulations where the N s sea quarks are in the p-regime, while the light valence quarks 
are in the e-regime, a setup that might be useful in mixed-action simulations. 

In either case A or B, we expect these predictions to match the ones of HMChPT in 
the limit — * oo, since this should recover a static limit of the valence quark. Indeed, the 
leading volume dependence in F^) an d Frg\ or Mm and Mm) can De shown to be associated 
to the light sector only and therefore should be independent of the heavy mass scale. We will 
explicitly show how this happens in section [71 

2.3. Pseudoscalar correlator 

Another interesting observable is the pseudoscalar density correlator as regards the finite 
volume dependence, because finite-size effects appear already at the leading-order as opposed 
to the correlator of the left current, where they appear first at NLO. 
The result at LO in the chiral expansion is the same for cases A and B: 



<3(0,0) 



Pi 



(2.33) 



y/v' 



C P (t) 



s 2 



P(t,M h ) {(U + ut) u ) + 2 



(2.34) 



2F 2 
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In this case, it is trivial to see that all the significant volume dependence comes from the 
zero- mode averages, which involve only the light sector. 



3. Static heavy-light mesons in finite volume HMChPT 

The effects of pion dynamics in the properties of static heavy-light mesons can be predicted 
in HMChPT |X6|, \T7\ 118] . Most calculations of chiral corrections have been done in infinite 
volume. The authors of [19] considered also chiral corrections in B parameters of neutral B 
meson mixing and heavy-light decay constants in a finite volume, but in the p-regime. We 
want to go further into the chiral limit by considering the e-regime for the light quarks. As 
far as we know, this regime has not yet been explored in HMChPT. However part of the 
technology we have used was developed in [20] to perform e-regime calculations in baryon 
ChPT. We present our NLO results for left-current correlators in the p-, e- and mixed regimes, 
and LO results for pseudoscalar density correlator. 

3.1. Formulation and conventions 

In the limit in which the mass of the heavy quark rrih goes to infinity it is expected that 
QCD simplifies. For example the interactions among the quark and the antiquark in a 
meson become spin independent, and if we consider processes in which only low momenta 
are involved, the heavy antiquark (or quark, if one prefers) can be decoupled. An effective 
field theory to analyse this situation can be built by rewriting the heavy-quark momentum 
p^ as: p^ = mfiVp + k^ and keeping only the leading term in the residual momentum k^/rrih- 
To recover the peculiarities of QCD, for example the chromomagnetic interactions, one has 
to rewrite the QCD Lagrangian as a series in powers of k^/rrih and consider also those terms 
that vanish in the rrih — > oo limit, up to the required degree of precision. 
In this work we have just considered the leading order in the above expansion. In such a 
case, the interactions with the pions are not able to modify the unitary velocity of the 
heavy-light mesons. We adopt a covariant representation, where the degenerate pseudoscalar 
and vector states are treated as a single field H which is usually labelled by v and the flavour 
a = 1, • • • , Ni of the light quark. In the Euclidean space we haveH 



3 While the formulation in Minkowsky space can be exhaustively found in the standard literature (see 
e.g. [21l [22]). we find useful to start from the beginning with the formulation in the Euclidean space. Notice 
however that for v/0 the Euclidean formulation is problematic [231 124) . and only the case v = (0, i) will be 
considered. 
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where P* and P represent respectively the vector and the pseudoscalar mesons, and P* 
satisfies: 

v ■ P* = . (3.3) 

The four- velocity v = (v, V4) satisfies the condition v 2 = — 1; the rest frame corresponds to 
v = (0, i). We use the conventional HQET normalisation of the states 

(H%\H b v ,) = 2v 4 (2ir) 3 5 vv ,6 ab , (3.4) 

according to which H fields have mass dimension -3/2. For simplicity, we drop the v label 
from here on. 

The Euclidean Dirac matrices are chosen to be Hermitean, 

tJ = 7am 75 = 75 = 71727374, (3.5) 
and satisfy the anticommuting relations 

{7am 74 = 2^- ( 3 - 6 ) 

The projector (1 — iv p ^ p )/2 in Eqs. (|3.1|3.2p retains only the particle component of the heavy 
quark. 

In a theory with N\ light quarks, and when dealing only with light mesons, one usually 
parametrises them with an SU(N[) matrix U = exp(2i£/F). If we rotate in flavour space the 
left (right) handed light quarks by a special unitary matrix L (R), the U field will transform 
like U — ► LUR) . As it is well known, when dealing with heavy- light mesons it is convenient 
to use the field yXJ to avoid that the parity transformation involves the pseudo-Goldstone 
boson field [El TH US] . y/U transforms like y/U -> L\/UW* or VU -» W\/UPJ, where W is 
a complicated function of R, L and the meson field £. Then H transforms as: 

H -> iJWt . (3.7) 

To write more easily a chiral invariant Lagrangian we build combinations of £, that like H, 
only transform with or W' under chiral rotations: 

V» = ^(Vu^d^x/U + VUd^), V p ^ WV^ + iWd^Wl (3.8) 

An = l -{Vu ] d^Vu -Vud^Vu ] ), a^wa p w ] . (3.9) 

Then, at leading order in 1/uih, a Lagrangian that is both Lorentz and chiral invariant is: 
^HMChPT = i^W a v^ ah + iV^H'] - ig n Tt[H a H b l5lu A b u a }. (3.10) 

The dynamics of the pseudo Nambu-Goldstone bosons is still given by the chiral Lagrangian 
inEq. (|2T2jh 
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From the kinetic part of Eq. (|3.10p one can extract the P and P* propagators. For v = (0, i) 
we obtain 

(P a (x)P b Hy)) = 5 ab V{x-y) (3.11) 
(P?(x)I**(y)) = 5 ab V(x-y)(6 flu -6^6 uA ), (3.12) 

where V(x — y) = ^5(x — y)9{x± — 7/4). See App. IA.1I for a more detailed discussion. 

The term of the Lagrangian in Eq. (13,101) proportional to g n represents the interaction of P, 

P* with an odd number of pseudo-Goldstone bosons. In particular, by expanding 

VU = e^ /F (3.13) 

we obtain the P*P^ and p*p*£ couplings 

Zhmchpt = - + ^du^ (P a] Pl* ~ PT ] P b ) + ^d v e a P^Pp*^v x (3.14) 

at leading order in the 1/rrih expansion. Note that the PPt; coupling vanishes because of 
parity. We adopt the convention 

£1234 = 1- (3.15) 

There are several determinations of g n on the lattice, in the quenched case |25[ [26| [27] and 
more recently in full QCD [2g [29l EQ] . 

A number of operators can appear at next-to-leading order in the chiral expansion |31j . 
however, if we omit contact terms, the only ones relevant to us are: 

5C H ] MChPT = -2viTr[HMH] - 2a' 1 Tr[HH]Tv[M}. (3.16) 

where A4 has been defined as: 

M = ^{VUM ] VU + VTJ ] MVU ] ). (3.17) 

The operator with the quantum numbers of the left current made of a heavy quark and a 
light antiquark with flavour index /, with the minimum power of H fields derivatives and 
mass insertions is: 

J^^Tr^^V] • (3-18) 

At leading order, the normalisation constant a is related to the pseudoscalar meson decay 
constant Fp and the corresponding mass Mp by the relation 

a = Fp^JlMp. (3.19) 

The vector meson decay constant is then given by 

F P , = MpFp, (3.20) 
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while for the masses one has Mp* = Mp. 

To represent the left current at NLO additional terms appear |31j : 



that absorb the UV divergences. 

In the static case, v = (0, i) the heavy-light left current correlator takes the form 



An 



(3.21) 



Q^cf\t) = J d 3 x (jl(x)J l u ($)), = (-V + 2(5^^4) , (3.22) 



where 



J 



-Tv^^P^UH) 1 ] . 



(3.23) 



Using this notation we isolate the time dependence in Cj (t) for later comparison with the 
mixed-regime result. We will use the index I = p to indicate the case where all light quarks 
are in the p-regime and I = e, where all are in the e-regime. Moreover, we will consider the 
case when some light quarks are in the p-, others are in the e-regime, and denote it by / = m. 
We are interested in the cases I = e and I = m, with e-regime valence quarks, which should 
match respectively the rrih — > oo limit of cases A and B in the ChPT computation. 

Similarly, at leading order in the momentum/mass expansion, the operator representing 
the pseudoscalar density is 



l5 H b (Vu bl + Vu 



t bl 



(3.24) 



where a is the normalisation factor defined in Eq. (|3.19p . In the case where all light quarks 
are in the e-regime, we will give the LO result for the correlator 



C P {t) = Jd 3 x(v l (x)V l (0) 



with 



P^Tr 
4 



75 H b (Vu m + Vu' 



t bl 



(3.25) 



(3.26) 



Note that we are using calligraphic characters to denote quantities calculated in HMChPT 
to distinguish them from the corresponding quantities of ChPT. 
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4. HMChPT in p-regime 

We consider HMChPT with Ni degenerate light quarks of mass m lying in the p-regime. 
Making use of the space integrals given in App. IA.2I we obtain, for t 7^ 0: 



Cf\t) = 0{t)%- exp (-AM^i) { 1 + 2m{ m + N lV3 ) 



+ 2F 2 L 2 { Nl Nj 



(P(t,M p )-P(0,M p ))(l + glj^ 



with M 2 = 2m£/F 2 and M p = \J M 2 + p 2 , while 



AM (p) = 2m(a 1 + Nia[) + 5; 



4F 2 L 3 V ' JVJ Y M 2 ' 



The function P has been already defined in Eqs. (|2.18|) . 



(4.1) 



(4.2) 



In dimensional regularisation ^2 p P(0, M p ) and ^ P -P(O) M p )M p 2 contain divergences, 
while ^ p M~ 2 is finite. To show this we rewrite: 



P(0,M p ) 



2M„ 



1 + 



e M p T _ X 



and define in s dimensions: 



G a>r (0,M) = — !— X)7T^ 

n*=i L * v (p + - 



M 2 ) T 



(4.3) 



(4.4) 



where in our case £1,2,3 = L,L^ = T. In the MS scheme we get: 

M 2 M 2 

G 4 ,i(0,M) = 2M 2 A(^) + -— 2 hi — + G^(0,M), 



(4vr) 2 /i 2 



G 3 3(0,M) = -8A(/i) 



1 

i^ 2 



M 2 
/^ 2 



ln^ + 1 ) + G^(0,M), 



G 3> i(0,M) 



In this expression A(/z) contains the divergence, 



A(/i) 



,4-d 



16vr 2 



1 1 

d-4 ~ 2 



(ln(4vr) + r'(l) + 1) 



(4.5) 
(4.6) 

(4.7) 
(4.8) 

(4.9) 
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while G\ r contains the finite volume dependence, which can be expressed as a series of Bessel 
functions: 



G&(0,M) ee _L^^_ l(M |z|) (4.10) 
G 3 y f (0,M) = ^L^^ (M|z|) (4.11) 



G^(0,M) = * Y, \ ^K_ h (M\z\) (4.12) 



where z = (niL%, • • • , n s L s ), {rii £ Z ; i = 1, • • • , s}. 
It can be shown that: 



-1^P(0,M P ) = G 4)1 (0,M), 
p 

1 v P(0,M p ) 1 rn , n | 1 y M p 3 r41 ^ 

L3Z^ M 2 " 2 G 3,|(°' M ) + X3 2. eMpT.T- ( 4 - 13 ) 

p P p 

Defining the renormalised coupling: 

m = r}\ r) +WiK») , (4-14) 
and requiring the cancellation of UV divergences we obtain, in agreement with [32^| 

fj + A^ 3 = §i( N i- ^) (1 + 3gl) . (4.15) 
Obviously one can also reproduce the infinite volume result by taking the limits T, L — > oo. 

5. HMChPT in e-regime 
5.1. Setup 

We consider now N[ light quarks lying in the e-regime. In this regime it is convenient to use 
the following parametrisation for the pseudo Nambu-Goldstone fields: 

U = U e^ (5.1) 

for which the integration measure is known up to NLO [33J and gives no contribution to our 
observables. Here £ contains the non zero modes of the pions and is a perturbative field £ ~ e. 



Note that there is no standard convention for the normalisation of the couplings 770 and 773. 
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The complication, in heavy-light mesons calculations, is that we need to express y/U as a 
function of \/Uo and £, up to e 2 corrections. The solution can be written in the form: 



U 



iA 



U ll + — 



(5.2) 



where A and B are Hermitian matrices {A is also traceless), respectively of order e and e 2 , 
linear and quadratic in the components of £. Imposing: 



we obtain the system of 2Nf — 1 equations: 



s/Tj) =jj + 0{e 



2\[U 



B + y/Uo'By/Ui 



4VU -AVUA + 2A 2 



(5.3) 



(5.4) 
(5.5) 



which can be solved in a particular system of coordinates for Uq. 

We have not found a simple way to solve the equations for general Ni , so we have considered 
the particular case of SU(2), that is iVj = 2. One convenient choice for this group is to use 
the hyper spherical coordinates: 



'Uq = cos ip + i sin ip sin 9 cos 4>a\ + i sin ip sin 9 sin <$>oi + i sin ip cos 9a^ , (5.6) 

where Oj are the Pauli matrices and the angle ranges are: 

ipe[0,Tv], 9£[0,tt], 0€[O,2tt]. (5.7) 

Note that to parametrise Uq we just need to extend the range of ip: £ [0, 2ir]. 

As usual, it is worth to perform the contractions of the non zero modes first and then 
perform the non perturbative integrations of ip, <p and 9 over the range specified by (|5.7p . 
The Haar integration measure to be used for the zero modes is, in hyper spherical coordinates: 



[VUo] 



7T- 



d A a5{a 2 



1) 



1 

2^ 



d^d9d(t> sin 2 2ip suit/ 



(5.8 



where a is defined through Uq = ao + i& • a. 



5.2. Left-current correlator 

If all light quarks are in the e-regime and for iVj = 2 we obtain 



cf(t)\ Nl=2 =9(t)^-exp(-AM^t 



1 + lj^2 L > T ) + ^(*» L ' T )) 



(5.9) 
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where 



H'(t,L,T) 



H{t,L,T) 




+ 73E p (*'IPl)- G (°> ) 



(5.10) 



and 



AM< £ ) 



(5.11) 



8F 2 L 3 ' 



The functions h± and P are defined in Eq. (|2,19p . (|2.18p . while the propagator G(0,0) is 
given in Eq. (EH3D . 

This expression contains no divergences in dimensional regularisation. It is interesting to 
stress the fact that the zero-mode integrals that contribute to various diagrams, nicely cancel 
in the sum of all contributions. In particular this means that the current correlator loses any 
dependence on quark masses close enough to the chiral limit, which also means no dependence 
on the topological sector. 

This result may be used to predict the behaviour of a correlator of left currents with the 
quantum numbers of the B meson, in a finite volume such that the u and d quarks are in the 
e-regime. 

In Fig. [H we show the ratio of the finite- volume to infinite volume correlator at t = 1 fm as 
a function of the volume for two boxes and two values of (g^ = and g n = 0.44, as recently 
computed on the lattice by [30]). Corrections are 0(3 — 4%) at 2 fm, and the dependence on 
£/„■ is mild. 

In Fig. [21 we show the time dependence of the correlator after factoring out the exp(— AM^ e 't) 
(we will see later in Sec. [8] that in any real fit to lattice data, AM' f ' would renormalise the 
static energy E stat ). 

5.3. Pseudoscalar density correlator 

For the pseudoscalar density, the result at the LO for arbitrary Ni is found to be: 



6. HMChPT in mixed-regime 
6.1. Setup 

In order to keep into account the effects due to the strange quark in heavy-light systems it 
is convenient to apply the power counting introduced in [8] and reviewed in section [2J To 



C l P (t) = ^9(t) <(Z7 + Ui) a ) + 2 



(5.12) 
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mi ~ e 4 , Nt = 2,N s = 




Figure 1: Ratio of C^ l (t = lfm) at fixed volume normalised to the oo volume result as a 
function of L for two boxes with T = L (solid) and T = 2L (dashed), and for g n = 0.44 (thick 
lines) [30] and g n = (thin lines). We have fixed S = (250 MeV) 3 , F = 90 MeV. 
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mi ~ e 4 , N t = 2, N s = 




l.l '- 



* ' \—i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i_ 

0.2 0.3 0.4 0.5 0.6 0.7 0.8 

t 

f 

Figure 2: AC^ 1 {t) / (a 2 exp(-AM( e )t)) as a function of t/T for T = L = 2 fm, and for 
9jt = 0.44 (thick line) [30] and g n = (thin line). 



18 



implement it in HQET, at least in the Ni = 2, N s = 1 specific case does not require more 
technology than the one introduced in the previous section. 

In practice all the steps described in the previous section must be applied again to the 
parametrisation given in Eq. (12. 7ft , 

6.2. Left-current correlator 

In this case the SU (3) vectorial symmetry is explicitly broken by the fact that 2 light quarks 
have mass m; lying in the e-regime while the one playing the role of the strange has a mass 
m s in the p-regime. This explains why the result is different depending on which light quark 
appears in the external line. As before we first consider the case in which I = 1, 2. This result 
represents the correlator of a left current with the quantum numbers of a B or a B* in the 
context of 2+1 light flavours. We obtained: 



This correlator will be matched with the predictions from the mixed ChPT, case B. 

In Fig. [3l we show the ratio of the finite-volume to infinite volume correlator as a function 
of the volume for two boxes and two values of g n in the mixed regime. We have set iVj = 2 and 
N s = 1. The corrections are qualitatively similar to those in the e-regime and quantitatively 
a bit larger. 

In Fig.HJ we show the time dependence of the correlator after factoring out the exp(— AM^ mi H). 

Another reason why HMChPT is useful is to predict the relation between observables re- 
lated to the B (B*) and the B s (B*). So we add for completeness also the results representing 
the correlator Cj^(t) of two left currents with the quantum numbers of a B s (or a B*). In 



.(m)l,2 

J 




(6.1) 



where M ria has been defined in Eq. (|2.30p and 




(6.2) 
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mi ~ e 4 , m s = 100 MeV, N { = 2, N s = 1 

1.10 
1.05 
1.00 
0.95 
0.90 

1 2 3 4 5 

L[fm] 

Figure 3: Ratio of Cj"^ 1 (t = 1/m) at fixed volume normalised to the oo volume result as 
a function of L for two boxes with T = L (solid) and T = 2L (dashed), and for g n = 0.44 
(thick lines) [30] and g n = (thin lines). The fact that the curves for T = 2L for g n = or 
0.44 nearly coincide is accidental. 
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mi ~ e 4 , m s = 100 MeV, N t = 2, N s = l 

i — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — r 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 

t 

T 

Figure 4: 4Cj Ti)1 (t)/(a 2 exp(-AM( mi ^)) as a function of t/T for T = L = 2 fm, and for 
g n = 0.44 (thick line) [30] and g n = (thin line). 
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this case we obtain 

Cf )3 (t) = #(t)^exp(-AAf( m 3) t ) |i + 2m s ( r?0 + r?3 ) 

iW(P(t,M ip )-P(0,M ip )) (l + 



where 



+ 



F 2 L 2 



alp 2 



p V V VbP 



AM^ = 2m s (a 1 +a[) + 



2F 2 L 3 



v M s 2 l v Ml 



(6.3) 



(6.4) 



Note that also in this case even though the various diagrams do depend on mi, the final result 
does not. 



7. Matching of HMChPT and ChPT 

Dominant finite-size effects in QCD are due to pion dynamics, since these are the lightest 
degrees of freedom. It is therefore expected that the finite-size scaling of heavy-light systems 
does not depend on the large energy scales related to the heavy quarks, ie. M^h or M^. 
This must be the case as long as those scales are significantly larger than L _1 . Whether 
these scales are much larger also than the QCD scale so that the static limit (HQET) is a 
good approximation, or not, should not matter a priori for the finite-size scaling properties, 
because the volume dependence arises from the propagation of the light degrees of freedom. 
The leading finite volume effects are therefore expected to come from the fact that the 
heavy meson can emit and absorb a pion. The probability for this to happen can however 
depend on the heavy mass scale. Close enough to the chiral limit, the masses of pseudoscalar 
mesons are suppressed by the spontaneous breaking of chiral symmetry, so, for example, we 
do not need to include the vector mesons in the effective theory, because they are much 
heavier and decouple. On the other hand, in the limit m/, — > oo pseudoscalar and vector 
mesons are degenerate, because the interaction between quark and antiquark inside the meson 
becomes spin independent, so they both need to be considered in HMChPT. The presence 
of heavy-light vector resonances can modify the finite volume effects indirectly by inducing 
unsuppressed contributions to pion/heavy- light meson scattering. We will see that indeed the 
finite-size corrections in HMChPT and mixed ChPT match up to corrections proportional to 

Real c and b quarks are somewhere in between these limits, where no effective description is 
very accurate. We may ask whether it is possible to give a description of finite size effects in 
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this intermediate regime. In particular, there might be other resonances to consider [34] . Us- 
ing general arguments it was shown in [35] that the pseudoscalar meson remains the lightest 
state for every value of the quark masses. Moreover if the heavy quark is in the non relativis- 
tic regime, we can say that the axial and scalar mesons (made of the same quark-antiquark 
couple) are heavier because they are in a higher angular momentum state. Experiments 
show that this peculiarity persists for heavy-light mesons whose heavy quark is a strange or 
a charm, the mass difference among the vector meson and the axial one being always of order 
of 400 MeV [36]. 

Finally the fact that exotic states may play a significant role is disfavoured by large N c argu- 
ments [37] saying that quark bilinears amplitude to produce them (like a qqqq) is suppressed. 
To sum up it seems plausible to consider a scenario in which the current correlator has two 
channels, a pseudoscalar and a vectorial one. The vectorial one could be integrated out for 
quark masses that are small compared to Aqcd- Indeed it is known [36] that while the 
Ks weight approximately 500 MeV the K*s weight approximately 900 MeV and while pions 
weight 140 MeV the ps weight 770 MeV. 

However the vectorial channel becomes more and more relevant as the mass of the heavy 
quark grows, because the mass difference between pseudoscalar and vector mesons dimin- 
ishes: if for K mesons it is about 400 MeV, for Ds it is about 150 MeV, while for Bs it is 
only 50 MeV. 

We consider now how the matching works in the two examples considered. Given any 
meson two-point function, the first point to realise is that a finite static limit is recovered 
after factorising out the leading e~ M ^, where M is the mass of the heavy meson and t is the 
temporal separation between the two mesonic sources. 

7.1. Pseudoscalar two-point function 

Let us start with the pseudoscalar correlator at LO, which is given in Eq. (|2,34p for the 
mixed ChPT case and in Eq. (15. 12ft for the HMChPT case. The first thing we observe is 
that the contribution of the zero modes, in particular the factor 



appear in both correlators. This shows that the zero modes contributions match in the two 
frameworks. 

Moreover, if we use the expansion 






-Mt 



(7.2) 



in Eq. (|2.34p we obtain 




S 2 



(7.3) 



AF 2 M h 
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After factorising out the exponential e~ Mflt we find that also the time-dependence matches 
exactly the one predicted by the HMChPT in Eq. (|5.12p . The matching of the coefficient 
gives 

**M k &r, ,7.4) 

2 F l M h m h 

which in the heavy quark mass limit ~ is consistent with the definition in Eq. (|3. 19p . 
a 2 /2 = FpMp. Since in the static limit there is no time dependence at LO, we have that the 
ratio of correlation functions at different volumes V\ and V2 is given by 

C l P (t)\vi _ ((Uo + ubu) + 2\vi , 7S) 
C l P {t)\v 2 {{U + Ul)u) + 2\ V2 

in both regimes of the heavy quark mass. 
7.2. Left-current two-point function 

We consider the matching for the Dirac components [i = v = 4 for which we have the ChPT 
result. 

What can be matched is the dependence of the correlators on the volume, that is L and 
T and the masses of the up, down and strange quarks, because these are explicit degrees 
of freedom in both effective theories. Moreover, since we only consider the static limit of 
HMChPT, we have to drop from the ChPT result those contributions that are suppressed by 
negative powers of m^. We expect that the e regime (J = e) HMChPT result should match 
to case (A) in the mixed-regime computation, while the I = m result in HMChPT should 
match case (B). 

• Case A: 

In order to match Eqs. (|5.9p and (|2.21j) . the L, T dependence must be the same in both 
cases. For the mixed ChPT framework, we split the contribution due to the heavy pions 
from the rest in Eqs. (|2.22p and (|2.23p and write: 



Ffa = F (m h , N t ) - 5 NVz - ^-J G(0,0) + O^ 1 ), (7.6) 



M( A) = M^mhNO + Oim^). (7.7) 

F and have absorbed the dependence on the heavy quark mass. The static limit 
of the mixed ChPT case in Eq. (|2.2ip is, for t > 0: 



Cf\t) F 2 M h 



exp(-M h t) 4 
where H (t, L, T) is the function of Eq. (|5.10|) 



1 + — \Ni- — \ Hit, L, T) 
2F 2 L 2 V Ni ' V ; 



(7.8) 
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For Ni = 2, the result is identical to the NLO prediction C$(t) in HMChPT (Eq. ([5T9]) ) 
with the following identifications: 



a = F P y/2Mp~ = F^2M h , g?T = 0. (7.9) 

The fact that at NLO we have to put g n = to match the two expressions reflects 
the fact that the vector meson is integrated out in the chiral theory. In HMChPT the 
vector and pseudoscalar are degenerate and therefore both are present. More generally 
we would expect that in the intermediate regime the finite size scaling of the current 
correlator behaves as 

C/ i/!|r ' 1 + ¥72 (Nt-^ltHfrLuTj+afrLuT^mh)) (7.10) 



Cf(t)\ v , 2 2F 2 L\ V N 



' f N-±- ) (H(t,L 2 ,T 2 ) + a(t,L 2 ,T 2 ,m h )) + ... (7.11) 



2F 2 L\ V Ni 
where V\ = L\T\ and V 2 = L^T 2 and 

lim mh _^ Qt(t, L, T, m h ) = 0, 
)im mh - too a(t,L,T,m h ) =&H'(t,L,T). (7.12) 

In the intermediate region the function a is unknown. However it should be possible 
to compute it including the leading 1/m/j corrections in HMChPT or even in ChPT 
including the vector resonance, as a function of the vector meson mass and coupling. 
We will not consider these regimes in the present work. Note however that for the value 
of g n obtained in a recent lattice computation in [30], g n = 0.44, the contribution of 
the term proportional to g\ (i.e. the difference between the thick and thin curves) in 
Figs. [T]|4]is not too large, and should decrease with decreasing m^. 

Case B: 



We have to follow the same steps as above, but in addition to L and T, we expect to 
reproduce also the dependence on m s , up to ni^ 1 contributions. In the mixed ChPT 
framework, we rewrite Eqs. (|2.3ip and (|2.32j) as: 

* = T 2, „ . _ lf„ 1 W_ N 
N 



F( B) = F (m h , N t + \N t - - G(0, 0) - -fG (0, M s ) + (7.13) 



* G{Q,M ns ) + %N s M 2 ss L A , 



2NNi 

o—2, \ 8M?„ M?N S (L 4 ~ 2L 6 ) 
Mf B) = M 2 h (m h ,N l + N s ) ss h p 2 ^. (7.14) 
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The mixed correlator in the static limit can then be written as: 



Cf\t) F 2 M h ( AM 2 ss M h N s (L A -2L 6/ 4 

■ exp t 



exp(-M h t) 4 F 2 

1 + + 2L 6 ) + ^ (jV, - j^j H(t, L, T) + 

+ ( Ns ^ (p (t ' Msp) " p (0 ' Msp)) + S ( p (*> M ^p) " P (°> M Vsp)) 

\ p 1 p / 

One can check straightforwardly that, for Ni = 2 and N s = 1, this coincides with the 
correlator Cj mi ^ 1,2 (i) computed in HMChPT (Eq. (16. ip ) with the identifications: 



a = FpyJ2Mp = Fy 2Mh, 

J = ^(4 } +24 p) ), 

o[ = -^(4 } -24 r) )- (7-15) 

Note that the above relations are among renormalised quantities. Apart from some 
finite volume effects due to the seap-regime quarks, which are exponentially suppressed, 
the volume dependence is identical to the one of Case A. So again we expect that for 
any value of m/,, Eq. (|7.1ip holds up to higher order chiral corrections and neglecting 
exponentially suppressed terms in exp(— M S L). 



8. Finite-size scaling of heavy-light mesons in lattice QCD 

As we have seen above the matching of finite-size effects of heavy-light correlators in HMChPT 
and ChPT works as expected. We are interested however in using these results to predict 
the finite-size scaling of these correlators computed in lattice QCD. On the lattice, we can 
include a relativistic or static heavy quark. In both cases we expect that for sufficiently large 
time separations: 

C l f(t) = EW*W°)>^ * C j(t) x ^exp(-Mt), (8.1) 

X 

where M is the lightest heavy-light meson mass M^i in the case of a relativistic heavy quark 

or the so-called static energy, E s t a t = Mhi — m h hi the lattice static limit. 

Note that the value of E s t a t is not predicted by HMChPT, however in general we can write: 

= E<£1 + AM (8.2) 
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where E stat is the value the static energy would have in the chiral limit, while AM contains 
the chiral corrections that are predicted by HMChPT, that we have presented for the various 
cases considered, in Eqs. ([33]) . (|5TTT]) . and (KT3|) . 

In practice this means that to fit a correlator evaluated with all the quarks in the e-regime 
using Eq. (|5.9|) one has to determine four parameters: a, F, E^ at and g w . It remains to 
be seen what the stability of such fits is in practice. The numerical challenge of extracting 
signals over the noise when computing propagators of heavy mesons is well known. Recent 
proposals to improve the situation have been discussed in [381 EH] • 



9. Conclusions 

We have considered the finite-size scaling of heavy-light mesons, composed of a light quark in 
the e-regime. We have computed the left-current and pseudoscalar two-point functions in two 
limiting regimes of the heavy quark mass: a small heavy quark mass such that the heavy-light 
meson can be treated in the mixed-regime of ChPT, and the static limit where HMChPT 
can be applied. We confirm the naive expectation that the dominant finite volume effects 
are induced by the emission/absorption of light pions, and are to a large extent insensitive 
to the value of the heavy quark mass. These results can be useful for matching lattice QCD 
and ChPT or HMChPT in finite volumes not sufficiently large compared with the Compton 
wavelength of the lighter pions. Our results can be used to consider also various partially- 
quenched situations. 



A. Space time integrations in HMChPT 

To obtain the charge correlators from the current ones, one has to integrate the current 
correlators over space. We report here the relevant results. 

A.l. The finite volume propagator in the rest frame 

In this section we want to obtain the propagator of the heavy-light mesons in Euclidean space 
and at finite volume. The propagator of HQET is obtained by writing the four-momentum 
of the heavy quark p^ as: p^ = ra^v^ + k^ and keeping only the leading term in the residual 
momentum k^. We consider here the rest frame, where v = (0,0, 0,i). For subtleties related 
to the Euclidean formulation for v / the reader can refer to \23\ [M] . In order to obtain 
the heavy quark propagator, one start from the Dirac quark propagator in coordinate space 
and we take the heavy quark limit, which is given by 

c M _ 1 f ,4 e^{-ip^ + m h ) h _ 1 + 74 [ + °° dp, e*4* _ 
boo[X) ~ (2k)* J P (pi + ml) ^ {t) ~ 2 2vr z(p 4 - im h ) ~ 

(A-l) 
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1 + 74 



6{t)e 



-m h t 



The projector (1 + 74)/2 retains only the particle content of the heavy quark, and for this 
reason the propagation in (|A-ip is forward in time. 

In the effective theory the exponential is factorised, and the static propagator at infinite 
volume is [2] 



Voo(x) 



1 



(2vr) 4 



~5(x)0(i) . 



(A-2) 



2i(p4 — ie) 

We now consider a finite box V = L 3 T with periodic boundary conditions. Analogously to 
(|A-ip . the finite- volume Dirac propagator in the heavy quark limit is given by 



S(x) -» S hq (t) 



1+74 1 



ip4,t 



that is, for < t < T, 



1 + 74 



1 



-m h T 



(A-3) 



(A-4) 



In the m/i — > oo limit, this reproduces the result of the infinite volume (|A-ip . Consequently, 
in the rest frame, the finite heavy volume propagator is 

1 



V(x) = -6(x)9(t), 



(A-5) 



which is exactly the propagator we obtain from the kinetic term of the HMChPT Lagrangian, 
Eq. (|3.10p . The heavy propagator has the same form as in infinite volume: this is not 
surprising, since it describes a static particle, which is not sensitive to the presence of a finite 
box. 



A. 2. Space integrals (p-regime) 



We present here the results for the integrals over space that are needed in finite-volume 
HMChPT when the light quark is in the p-regime. V{x) represents the static propagator, 
Eq. ()A-5j) . while G(x,M) is the pion propagator of Eq. (j2.14[) . The function P(t,M) is 
defined in Eq. (|2TT8|) . 



A 2 (t,M) 
A 3 . a (t,M) 



d i xV{x) = h{t)- 



d 3 xV(x)G(x,M) = 



2L 3 



^P(t,M p 



d 3 xd 4 zV(x - z)V{z)d Xa G{x - z,M) 



6{t) 



^P(t,M p )-G(0,M) 



L 3 



(A-6) 
(A-7) 



(A-8) 



2S 



d 3 xd 4 zd 4 wV(x - z)V(z - w)V{w)d Za d Wl3 G{z -w,M); (A-9) 



if a / /? 
1 



A 4;a/3 (t,M) 
A 4 . a(3 (t,M) 

A 4 . M (t,M) 

a 

A. 3. Space integrals (e-regime) 



G(0,M)--i^P(t,M p ) 



M2 



+ P(t,M p )-P(0,M p 



8L 3 ' 

d 3 xd 4 zy(x - z)V(z) = ^tO(t); 
d 3 *d A zd A wV{x - z)V(z - w)V{w) 
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t 2 8(t). 



(A-10) 
(A-ll) 

(A-12) 

(A-13) 
(A-14) 



In the e-regime the integrals to be computed are the same as above, with G(x, M) substituted 
by G(x, 0) defined in Eq. (|2.15p . We will denote the corresponding integrals by A n (t) instead 
of A n (t,M). We have obtained: 



Mt) 

A 3a {t) 
A 4 . 44 (t) 

^ A 4 - a a{t) 

a 



2L 3 



Ja4- 



6(t) 



pfcif^ + pE^'ipi)-^ ' ) 



T 
Z 3 ' 



1 t 2 

—e(t)-; 

16 1 V 



(A-15) 
(A-16) 

(A-17) 

(A-18) 
(A-19) 



if a ^ /? . 

The function h\{t/T) is defined in Eq. (12. 19ft . Notice that no new integrals have to be 
considered in the mixed-regime case. 
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